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Abstract. We recently put forth a new fundamental lattice Hamiltonian based on an 
underlying picture of electrons and deuterons as elementary Dirac particles. Within 
this model there appears a term in which lattice vibrations are coupled to internal 
nuclear transitions. This is interesting as it has the potential to provide a connection 
between experiment and models that describe coherent energy transfer between two- 
level systems and an oscillator. In this work we describe a calculation of the coupling 
matrix element in the case of the deuteron based on the old empirical Hamada- Johnston 
model for the nucleon-nucleon interaction. The triplet S and D states of the the 
deuteron in the rest frame couples to a singlet P state through this new interaction. 
The singlet P state in this calculation is a virtual state with an energy of 125 MeV, 
and a coupling matrix element for z-directed motion given by 2.98 x 10 -3 MjcP z . 



1. Introduction 

We recently obtained a new Hamiltonian for a lattice which includes interactions with 
internal nuclear degrees of freedom pQ . We started with a Dirac model for electrons and 
for nucleons on equal footing, then allowed the nuclei to be described by a finite basis 
expansion, and finally developed an approximation appropriate for low nuclear velocity. 
Our original motivation for this was to obtain a model capable of describing the mass 
shift associated with excited nuclear states in a configuration interaction calculation. 
However, the new model unexpectedly contains a new coupling term which provides for 
a direct interaction between lattice vibrations and nuclear transitions. 

On the face of it, this new interaction term looks like it should allow for the strongly- 
coupled transitions that we have sought [2] in connection with the generalized lossy spin- 



Coupling between a deuteron and a lattice 2 

boson models [31 HJ El El El IE] that we proposed to account for some of the anomalies 
(such as excess heat P QUI HH [12] and collimated x-ray emission [13], [HI [151 EH ITT] ) 
that have been reported in experiments in recent years. What is needed at this point 
is an explicit calculation for some nuclear system to see how it works, what states are 
coupled to, and how big the interaction is. 

In general nuclear structure models are much more complicated than atomic 
structure models due to the more complicated strong force interaction. We would like 
to work with empirical strong force models derived from scattering experiments and 
few-body bound state binding energies. In recent years these models have achieved 
impressive results [HI [19]; however, some of these strong force models involve a fair 
amount of work to implement. If we go back a few decades we can find simpler versions 
of strong force models that are easier to work with, and are sufficiently accurate to 
clarify the issues of interest here. In the computations that follow we will focus on the 
old Hamada- Johnston potential model [20]. Without question the simplest compound 
nucleus which should show the effects of interest is the deuteron, and so we will focus 
on this system in what follows. 

In this formulation we have modeled the nucleons as elementary Dirac particles. 
As nucleons are made up of strongly interacting quarks, we know that they are not 
elementary Dirac particles. To do better in the case of coupling with the deuteron, we 
would require a description in terms of the six constituent quarks. We would expect 
from such a model a coupling matrix element likely somewhat different from what we 
calculated in this work. Even so, it makes sense here to pursue this simpler deuteron 
model based on simple Dirac nucleons as a step forward in the modeling process. 



Coupling between a deuteron and a lattice 
2. Basic model 



3 



In a recent paper we discussed the derivation of a finite basis approximation for a moving 
nucleus in the many-particle Dirac model which leads to the new coupling that we are 
interested in. We begin with the (relativistic) finite basis model that we obtained. 

2.1. Finite basis approximation 

In [1] we developed finite basis eigenvalue relations in the form 



E Cn = V /(M n C 2 ) 2 + C 2 |P| 2 C n + H nmCm (1) 

where the off-diagonal matrix elements were written as 

H nm — CX-nm ' (cP) 4" V nrn (2) 

Here V nm is the relative coupling matrix element 



O (3) 



We defined a nm as 

/ rn ■ \ 

$m> (4) 



Otnm. = <$ 



n 



x ^ rrij 



The notation for the two-body version of the problem is a bit different than what 
we used for the many-particle problem. It is useful to recast the relative matrix element 

as 



V — < $ 



$m) (5) 



(a 2 - cti) ■ cp + fiirrixc 2 + /3 2 m2C 2 + V(r) 

The relative part of the off-diagonal matrix element corresponds to the rest frame 
interaction terms, which might come about from strong force interactions as in the 
development above for the nonrelativistic deuteron problem. What is new is the coupling 
with the center of mass momentum P that appears in a nm ■ (cP). We are interested in 
these new matrix elements. 

2.2. Nonrelativistic reduction of the new interaction term 

In pp we discussed the reduction of the new interaction matrix element to the 
nonrelativistic case. The results can be expressed as 
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As above, this is written for the many-particle problem, and we wish to recast it in 
terms of the two-body problem; we may write 



^ M 3 

3 



<I>, > -)• 
1 



(£ - Mc 2 
2Mc 2 













"-) 


(*/ 


p-P 


\m 2 


mi/ 





2Mc 2 

+(*/ 
1 



>i • cP) 



1 



2mic : 



■V(r) 



(o-i • cp) 



^i • cpj 



2m 1 c 1 



■V(r) 



:<ti • cp; 



2Mc 2 



,0-2 • CP) 



>2 • CP) 

1 V(r) 



V(v) 



2m 2 c 2 



2m2C 2 

0- 2 • CP) 



>2 • Cp) 



(7) 



^.c?. Equal mass approximation 



It is possible to split up this new interaction term into a contribution that takes the 
nucleon masses to be equal, and a small correction term that depends on the difference 
between the nucleon masses. In what follows our focus will be on the larger equal mass 
terms, which is equivalent to making an equal mass approximation. In this case we may 
write 



\- rrij „ 
^ M 3 



+ 



1 



°i ■ cpj 
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where we have assumed that 



V{r) 
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mi = m 2 = m a 



(8) 



(9) 



In this approximation there is no longer an explicit dependence on the state energy E. 
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2.4- Nonrelativistic approximation 

It is possible to develop a nonrelativistic approximation using 



^(M n c 2 ) 2 + c 2 |P| 2 -> M n c 2 + 



2M„ 



(10) 



In this finite basis model that includes center of mass dynamics to lowest order 

could be developed starting from a nonrelativistic Hamiltonian of the form 



H = Mc 2 + M 7 + ^ + y(r) + a-cP 



2M 2/i 

where M is an operator that returns the rest mass energy of the nuclear state, and 
where a • cP in the equal mass approximation is 



a cP ->■ 



2Mc 2 / \2m av c 2 



{(T 2 -Cp)V((T 2 -Cp) + ((T r Cp)V((T r Cp) 



- (<ri • cP)V (a, ■ cp) - (eri • cp)f( - 1 • cP) 



(12) 
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We know from the literature that the deuteron at rest can be modeled using a triplet S 
and triplet D state, since the tensor interaction mixes the two. Since the kinetic energy 
and potential terms preserve J and Mj, each of the triplet S states mixes with a triplet 
D state that has the same J and Mj. 

3.1. Mixing with 1 P states 

The new interaction term causes these states to mix with singlet P states. In general, the 
new term does not preserve Mj, so that we would require a finite basis approximation 
that distinguishes the different sublevels. However, it is possible to focus on a special 
case of the new interaction which does preserve Mj. This occurs if we restrict our 
attention to 

P = i z P z (13) 

We find in this case that mixing occurs for \Mj\ = 1, but not for Mj = 0. In response, 
we might write 

\p = *3 S + *3 D + ^1/lp (14) 

with the understanding that 

^ lp (J M = 0) ->• (15) 

3.2. Basis state construction 

Nuclear state construction is usually carried out in the isospin scheme, with 
antisymmetry enforced through the application of the generalized Pauli principle. The 
two-body problem is particularly simple in this regard, with spin, isospin and spatial 
components restricted to being either symmetric (s) or antisymmetric (a); we may write 



for the three states 

*3 S = R{s)S{s)T{a) (16) 

^ 3D = R(s)S(s)T(a) (17) 

tfip = R(a)S(a)T(a) (18) 



The antisymmetric spin and isospin terms S(a) and T(a) are singlets, and the symmetric 
spin and isospin terms S(s) and T(s) are triplets. 
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The S state is a triplet spin state, so we may write it as 
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#3 S = #3 S (S=l,Ms;T = 0,Af r = 0;Z = 0,m = 0) 
= ^V 00 (M) \l,M s ) s |0,0) r 



(19) 



The Yi m (0, <p) are spherical harmonics; we choose / = and m = since we are working 
with an S state. The \S, M$)s are spin functions for the neutron and proton spins; the 
|T, M T ) T are isospin functions, and we have used an isospin singlet function here. 

3.4- Triplet D state 

We can develop a D state by applying the tensor ,512 operator on an S state. This 
approach was used early on as a convenient way of generating few-body wavefunctions 
for variational calculations in nuclear physics. We may write 



*3 D 
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This construction is convenient since 

v(r) 
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2^ 



3D 



3.5. Singlet P state 

The singlet P state for a particular calculation can be specified using 

tfip = tfip(S = 0,M s = 0;T = 0,Mr = 0;Z = l,m) 
w(r) . 



-Yi r 



\0,0)s |0,0) r 



Including an i here leads to real coupling coefficients in what follows. 



3.6. Normalization 

We can evaluate the normalization integral for these states simply; we write 



(20) 



(21) 



(22) 



= (^3 S |*3 S ) + (*3 D |*3 D ) + (*lp|*lp) 

roo 

/ \u(r)\ 2 + \v(r)\ 2 + \w(r)\ 2 dr 
Jo 



(23) 
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3. 7. Expectation value of Hamiltonian terms 

We are interested in developing coupled channel equations that include the new 
interaction. For the problem in the rest frame, this is most easily accomplished by 
developing an expression for the total energy and then using the variational principle. 
We can use the same basic approach here for the moving frame version of the problem. 
We begin with 



! + + CP), 
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+ v 
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(a-cP) 2 ^3 D ) + (^ S (a-cP) 2 ^i P ) + (^3 D (a-cP) 2 *i P ) (24) 



3.8. Diagonal matrix elements 

We can evaluate the diagonal matrix elements directly using Mathematica to obtain 
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(25) 



3^(r)+64 t (r)-3< 5 (r)-3<(r) 



v(r)dr(26) 



2(i dr 2 



^v° c s {r)-2v° L s L {r) 
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3.9. Off- diagonal potential matrix elements 

In the case of the Hamada- Johnston potential, there occur off-diagonal matrix elements 
between the triplet S and singlet D states, which are given by 



V 



V 



^3 S ) = - 3y/8 J™ u(r)v$v(r)d,7 



(28) 



The superscript et in the associated potentials here is connected with the even 
triplet channel, since the Hamada- Johnston potentials are fit for the different channels 
separately. 
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3.10. Off-diagonal matrix elements for the new interaction 

For the off-diagonal matrix elements of the new interaction, we have used Mathematica 
to compute 



<$ip|(a-cP) z |$ 3s > = Mj 
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(31) 
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We have specified a finite basis problem with three channels, which would produce 
three complicated coupled-channeled equations if we decided to treat the different basis 
states on equal footing. However, since the momentum P that we are interested in for 
applications of this model is small, the triplet S and D channels are then best considered 
to constitute the unperturbed deuteron problem, and the singlet P channel will contain 
the weak response of the deuteron to the a • cP perturbation. 

In this case, it seems appropriate to develop the coupled triplet S and D channels 
consistent with the rest frame deuteron problem. Once the associated wavefunctions are 
known, then we can use them to approximate the occupation of the singlet P channel. 

4-1. Rarita-Schwinger equations 

Given the approach outlined above, we can optimize the channel wavefunctions u(r) 
and v(r) by minimizing the rest frame energy 
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The minimization of this rest frame energy leads to the constraints 
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(35) 



where E r is the relative energy. We recognize these as the Rarita-Schwinger equations 
based on the Hamada- Johnston potential model. 
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Figure 1. Numerical solutions of the Rarita-Schwinger equations for the deuteron 
using the Hamada- Johnston potential. The solution for the S state \u(r)\ is shown in 
blue; the solution for the D state [v(r)] is shown in red. 



4-2. Rest frame triplet S and D channel wavef unctions 

We have solved the Rarita-Schwinger equations to obtain the channel wavefunctions 
plotted in Figure [H The triplet S channel wavefunction u{r) is larger and extends out 
to a relatively large radial separation, and the triplet D channel wavefunction v(r) is 
smaller and localized to much smaller radial separation. We can see the effect of the 
hard core potential in the zero boundary condition at the cut off radius. 

4-3. Optimization of the singlet P channel 

In the perturbation theory approach outlined above, we can approximate the occupation 
of the singlet P channel in terms of known triple S and D channel wavefunctions. The 
associated constraint on the channel wavefunction can be written as 
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Figure 2. Numerical solution for the normalized singlet P radial wavefunction. 
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We have solved this equation numerically assuming that u{r) and v(r) are fixed solutions 
of the Rarita-Schwinger equations, and the resulting normalized solution for w(r) is 
shown in Figure [2j 



4-4- Equivalent two-level model parameters 

From the computation outlined above we can derive an equivalent two-level system 
model in the form 



E 



Cl 



c 2 



H21 H22 



We compute 



H22 



- 2.245 MeV 
125.4 MeV 



C-2 



(37) 

(38) 
(39) 



Coupling between a deuteron and a lattice 



14 



H 12 = 



21 — 



2.98 x 1(T 3 Mj(cP z ) 



(40) 



The off-diagonal coupling matrix elements are somewhat smaller than we were hoping 
for, and future work will be needed to understand if this coupling is sufficiently large 
to account for experimental results. In addition, we have found that these off-diagonal 
matrix element depend on the nuclear spin, which suggests that the system may respond 
to net spin alignment. 
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We recently proposed a new fundamental Hamiltonian for condensed matter lattice 
problems that includes coupling to nuclear internal degrees of freedom. From our 
perspective this new coupling seems to be what is needed to account for the excess 
heat effect in the Fleischmann-Pons experiment. What has been needed in order to 
evaluate the models that result is an estimate for the coupling matrix element. 

The development of an estimate for this matrix element is challenging for a variety 
of reasons. We have presumed in the derivation of the fundamental Hamiltonian that 
it is sufficient to model the nucleons as elementary Dirac particles. However, we know 
that nucleons are composite particles made up of quarks and gluons, and that it is 
unlikely that using a Dirac model as we have done is going to give accurate results. To 
do better we probably need to go back and develop a better fundamental Hamiltonian 
based on quarks and electrons. If it is possible to obtain reasonable nucleon models from 
empirical potentials, then we may be able to develop a better estimate for the deuteron 
coupling matrix element. Working directly with bound state QCD at this stage does 
not seem to be an attractive option. 

Once we have decided on the simpler model that adopts an elementary Dirac particle 
model for nucleons, then it is an issue of whether to use a relativistic or nonrelativistic 
model, and further it is an issue of what potential to use. Since these computations 
involve a fair amount of work, it seemed sensible to adopt a nonrelativistic model since 
it is simpler, and to work with an older relatively simple nuclear model. The Hamada- 
Johnston potential fits the bill in this regard, as it is sufficiently simple that we are 
able to complete a calculation in relatively short order. Perhaps the most work in this 
computation was the evaluation of the spin, isospin, and angular momentum algebra; 
for this we relied on brute force Mathematica calculations. 

In the end, we have developed a model for the coupling between the different nuclear 
spin states of the ground state deuteron and lattice-induced coupling to a highly-excited 
singlet P virtual state. The energy of this virtual state is about 125 MeV in this model, 
which is consistent with our expectations. The coupling matrix element fell short of 
what we had hoped for by about an order of magnitude. We will need to clarify in 
future calculations if this is sufficiently large to be relevant to experimental results. 

The coupling matrix element in this model is proportional to Mj, which is 
interesting in connection with the reported dependence of excess heat on the strength of 
an applied magnetic field. Since the matrix element is proportional to Mj, there is the 
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potential for a larger coupling if the deuteron spins can be aligned. We are interested 
in pursuing this possibility in future work. 
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